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I. INTRODUCTION 

The study of the production properties of the 
Z boson began in 1983 with its discovery by the 
UA1 and UA2 collaborations at the CERN pp col- 
lider jl||- Together with the discovery of the W 
boson [3|4j earlier that year, the observation of the 
Z boson provided a direct confirmation of the unified 
model of the weak and electromagnetic interactions, 
which, together with QCD, is now called the stan- 
dard model. Since its discovery, many of the intrin- 
sic properties of the Z boson have been examined 
in great detail via e + e~ collisions at the LEP e + e _ 
collider at CERN || . The mass of the Z boson mea- 
sured at LEP and the SLC e + e~ collider at SLAC, 
known to better than one part in 10 4 [||, is one of 
the most precisely measured parameters in particle 
physics. 

LEP experiments have focused on the intrinsic 
properties of the Z boson, examining the elec- 
troweak character of its production and decay in 
e + e~ collisions. At the Tevatron, where the Z boson 
is produced in pp collisions, its production proper- 
ties are presumably characterized by QCD. Since the 
electroweak properties of the Z boson are not cor- 
related with the strong properties of its production, 
the Z boson can therefore serve as a clean probe of 
the strong interaction. Also, the large mass of the 
Z boson assures a large energy scale (Q 2 w Mf ) 
for probing perturbative QCD with good reliability. 
The measurement of the cross section as a function 
of transverse momentum (da/dpx) of the Z boson 
provides a sensitive test of QCD at high-Q 2 . In this 
article, we describe a measurement of da/dpr of the 
Z boson using the e + e~ decays of the Z M. 

In the parton model, at lowest order, Z bosons 
are produced in head-on collisions of qq constituents 
of the proton and antiproton, and cannot have any 
transverse momentum. Consequently, the fact that 
observed Z bosons have finite pr is attributed to 
gluon radiation from the colliding partons prior to 
their annihilation into the Z boson. Gluon radia- 
tion within the color field of the proton or antipro- 
ton increases in proportion to the time available for 
such annihilation, which is proportional to the in- 
verse of the energy scale for the process (l/Q) [||. 
The radiated gluons carry away transverse momen- 
tum from the annihilating quarks and momentum 
conservation requires that this be observed in the 
Pt of the Z boson. Thus, one expects that the ob- 
served transverse momentum distribution of any di- 
electron system (produced at a scale Q ss M ee ) will 
broaden as a function of Q. This is, indeed, the ef- 
fect observed. At M ee « 10 GeV, the typical pr 
for Drell-Yan pairs § is about 1 GeV gcj. For W 
boson production (Q sa 80 GeV), the average pr is 
about 5 GeV H]. For Z boson production (Q as 91 
GeV), the average pt is about 6 GeV pjfl. 



In general, the differential cross section for pro- 
ducing the state V is given by: 



dp'Ldy 



dx l dxjf{x i )f{xj)- 



dp 2 T dy ,^ 



where pr and y are the transverse momentum and 
the rapidity of the state V; Xi and Xj are the mo- 
mentum fractions of the colliding partons; f{xi) and 
f(xj) are the parton distribution functions (pdf's) 
for the incoming partons; and <7y_>y is the partonic 
cross section for production of the state V, in our 
case, the Z boson. The subscripts i and j denote 
the contributing parton flavors (i.e., up, down, etc.) 
and the sum is over all such flavors. 

In standard perturbative QCD (pQCD), one cal- 
culates the partonic cross section by expanding in 
powers of the strong coupling constant, a s . This 
procedure works well when p T ~ Q 2 . However, as 
Pt — > 0, correction terms that are proportional to 
a s ln(Q 2 /p T ) become significant for all values of a s , 
and the cross section diverges at small px- Phys- 
ically, the failure of the calculation is due to the 
presence of collinear and \ow-px gluons that are not 
properly accounted for in the standard perturba- 
tive expansion. This difficulty is surmounted by re- 
ordering the perturbative series through a technique 
called resummation ||JRJ-|1|. 

In its final form, the differential cross section is 
calculated as a Fourier transform in impact param- 
eter, b, space: 



d 2 a. 



dp\dy 



d 2 b e l P T b W(b,Q) + Y(b, Q) 



(2) 



where W(b, Q) contains the results of resumming the 
perturbative series, and Y(b, Q) adds back to the 
calculation the pieces that are perturbative in a s , 
but are not singular at pt — ||. 

Although the resummation technique extends the 
applicability of pQCD to lower values of px, a more 
fundamental barrier is encountered when px ap- 
proaches Aqcd, and pQCD is expected to fail in 
general. In this region, we expect non-perturbative 
aspects of the strong force to dominate the produc- 
tion of the vector boson. This implies that W(b, Q 2 ) 
in Eq. ^| is undefined above some value of b = b max . 
To extend the calculation to px = 0, the following 
substitution is made: 



W(b,Q) -» W{K 1 Q)e- SNP ^ b ' Q) 



(3) 



where 6* = bj + (6/6 ma x) 2 - This effectively cuts 
off the contribution of W(b, Q) near 6 max , leaving 
the differential cross section dominated by Snp, 
where SNp(b,Q) is called the non-perturbative Su- 
dakov form factor. Snp has the generic renormal- 
ization group invariant form S 
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S N p(b,Q) = 

fci(M<) + hx{b, Xj ) + h 2 (b) ■ H^-J ( 4 ) 

where X4 and Xj are the momentum fractions of the 
annihilating quarks; Q Q is an arbitrary momentum 
scale; and h\{b 1 x), h^ib) are phenomcnological func- 



tions to be determined from experiment fll5|, 17| . |l8j . 
The fact that hi (b) lacks any dependence on the mo- 
mentum fractions of the incoming partons has led 
to speculation that it may contain some deeper rel- 
evance to the gluonic structure of the proton [^o| . 

The current understanding of the pt distribution 
for Z bosons uses fixed-order perturbative calcu- 
lations [leading-order (LO) or next-to-leading-order 
(NLO)] to describe the high-py region, and resum- 
mation calculations of the perturbative solution to 
describe the low-pr region. The resummation cal- 
culation fails at large-pT (of the order 50 GeV) due 
to large terms missing from the calculation resulting 
from the px — * approximation. An ad hoc "match- 
ing" criteria is invoked to decide when to switch from 
the resummed calculation to the fixed-order calcu- 
lation, which is considered to be robust at large-p-r- 
Additionally, a parameterization of Eq. |] is invoked 
to account for non-perturbative effects at the lowest 
Pt values which are not calculable in perturbative 
QCD. 

In our measurement of the px distribution, we re- 
strict the invariant mass of the dielectron system to 
be approximately the mass of the Z boson, where 
the Z resonance greatly dominates dielectron pro- 
duction. The remaining contribution is due almost 
entirely to production of e + e _ pairs via the photon 
propagator (Drell-Yan process), which is coherent 
and interferes quantum mechanically with Z boson 
production. Other processes also contribute to in- 
clusive dielectron production in pp collisions, e.g., ti 
and diboson production, however, these are incoher- 
ent with Z boson production and their overall rate 
is negligibly small. 

Besides being of intrinsic interest in the study 
of QCD, precise understanding of Z boson produc- 
tion in pp collisions has important practical benefits 
for other measurements with electrons in the final 
state. The phenomenology used to describe Z bo- 
son production is applicable to W, Z, and essentially 
all Drell-Yan type processes. In the low-pr region, 
where the cross section is highest, uncertainties in 
the phenomenology of vector boson production have 
contributed to the uncertainty in the measurement 
of the mass of the W boson (M w ) Addi- 
tionally, diboson, top quark, and Higgs boson pro- 
duction all have single and dielectron backgrounds 
from W and Z boson production that will be more 
constrained through a precise measurement of Z bo- 
son production properties. 

Despite larger statistical uncertainties relative to 
W boson production (there are ~10 times more 



W — > ev than Z — > ee events produced at s/s =1.8 
TeV), the Z boson provides a better laboratory for 
evaluating the phenomenology of vector boson pro- 
duction. The measurement of the transverse mo- 
mentum of the e + e~ pair (pf?) does not suffer from 
the same level of experimental imprecision as the 
measurement of p T v because the latter relies on the 
determination of the total missing transverse mo- 
mentum in the detector ($ T ), which has inherently 
higher systematic uncertainties. The typical resolu- 
tion in Pt is about 1.5 GeV compared with 4-5 GeV 
for pjf, and the p^ resolution is approximately flat 
as increases, whereas it continues to degrade for 

Previous measurements of the differential cross 
section for Z — > ee production in pp collisions have 
been limited primarily by statistics. The UA2 f23|| 
collaboration analyzed 162 events, concluding that 
there was basic agreement with QCD, but that more 
statistics were needed. In the 1988-89 run at the 
Tevatron, CDF jl2| analyzed 235 dielectron events 
and 103 dimuon events, making similar conclusions. 
Our study is based on a total of about 6400 events. 
We determine the px distribution for the Z boson 
and use our results to constrain the non-perturbative 
Sudakov form factor. We then remove the effects of 
detector smearing and obtain a normalized differen- 
tial cross section da/dpr- 

We present a brief description of the DO detector 
in the next section. We then present the selection 
procedure for our data sample. The selection effi- 
ciency (Section ll\I) , kinematic and fiducial accep- 
tances (Section |V| ), contributing backgrounds (Sec- 
tion [Vl), fit for non-perturbative parameters (Sec- 
tion |VTIID , and the smearing correction (Section [Tx] ) 
are all discussed in turn. These individual compo- 
nents are combined (Section |x|) to obtain the final 
differential cross section, which is compared to pre- 
dictions from QCD. 



II. EXPERIMENTAL SETUP 

The DO detector consists of three major subsys- 
tems: a central detector, a calorimeter (Fig. [j]), and 
a muo n sp ectrometer. It is discussed in detail else- 
where [g4| . We describe below only the features that 
are most relevant for this measurement. 



A. Conventions 

We use a right-handed Cartesian coordinate sys- 
tem with the z-axis defined by the direction of the 
proton beam, the a;-axis pointing radially out of the 
Tevatron ring, and the y-axis pointing up. A vector 
p is then defined in terms of its projections on these 
three axes, p x , p y , p z . Since protons and antiprotons 
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in the Tevatron arc unpolarized, all physical pro- 
cesses are invariant with respect to rotations around 
the beam direction. It is therefore convenient to use 
a cylindrical coordinate system, in which the same 
vector is given by the magnitude of its component 
transverse to the beam direction, pr, its azimuth 0, 
and p z . In pp collisions the center of mass frame 
of the parton-parton collisions is approximately at 
rest in the plane transverse to the beam direction, 
but has an unknown boost along the beam direc- 
tion due to the dispersion of the parton momentum 
fraction within the interacting proton and antipro- 
ton. Consequently, the total transverse momentum 
vector in any event (Ifix) must be close to zero, and 
can be used to reject background from events that 
have neutrinos in the final state. We also use spher- 
ical coordinates by replacing p z with the colatitude 
9 or the pscudorapidity 77 = — In tan (0/2). The ori- 
gin of the coordinate system is, in general, defined 
as the reconstructed position of the pp interaction 
for describing the interaction, and the geometrical 
center of the detector when describing the detector. 
For convenience, we use natural units (H = c = 1) 
throughout this paper. Additionally, we use "pt" to 
refer to the transverse momentum of the Z boson or 
objects which mimic the Z, e.g., background events 
in which the momenta for the objects considered to 
form the iake-Z are added together to generate a pr 
value. Deviations will be noted with an appropriate 
superscript. 



D0 LIQUID ARGON CALORIMETER 




FIG. 1. A cutaway view of the D0 calorimeter and 
tracking system. 



B. Central Detector 

The central detector is designed to measure the 
trajectories of charged particles. It consists of a ver- 



tex drift chamber, a transition radiation detector, a 
central drift chamber (CDC), and two forward drift 
chambers (FDC). There is no central magnetic field, 
and D0 therefore cannot distinguish particles by 
their electric charge, with the exception of muons 
which penetrate the outer toroidal magnets. Con- 
sequently, in the rest of this paper, the term elec- 
tron will refer to either an electron or a positron. 
The CDC covers the detector pseudorapidity region 
|?7dct| < 1-0. It is a drift chamber with delay lines 
that give the hit coordinates along the beam direc- 
tion (z) and transverse to the beam (r, <f>). The FDC 
covers the region 1.4 < |r?dct| < 3.0. 

C. Calorimeter 

The sampling calorimetry is contained in three 
cryostats, each primarily using uranium absorber 
plates and liquid argon as the active medium. There 
is a central calorimeter (CC) and two end calorime- 
ters (EC). Each is segmented into electromagnetic 
(EM) sections, a fine hadronic (FH) section, and 
coarse hadronic (CH) sections, with increasingly 
coarser sampling. The entire calorimeter is divided 
into about 5000 pseudo-projective towers, each cov- 
ering 0.1x0.1 in 77 x (p. The EM section is seg- 
mented into four layers which are 2, 2, 7, and 10 
radiation lengths in depth respectively. The third 
layer, in which electromagnetic showers reach their 
maximum energy deposition, is further segmented 
into cells covering 0.05x0.05 in 77 x <p. The hadronic 
sections are segmented into four (CC) or five (EC) 
layers. The entire calorimeter is 7-9 nuclear interac- 
tion lengths thick. There are no projective cracks in 
the calorimeter, and it provides hermetic and nearly 
uniform coverage for particles with |?7det| < 4. 

D. Trigger 

Readout of the detector is controlled by a multi- 
level trigger system. The lowest level hardware trig- 
ger consists of two arrays of scintillator hodoscopes, 
which register hits with a 220 ps time resolution and 
are mounted in front of the EC cryostats. Particles 
from the breakup of the proton and the antiproton 
produce hits in hodoscopes at opposite ends of the 
CC, each of which are tightly clustered in time. At 
the lowest trigger level, the detector has a 98.6% 
acceptance for W/Z boson production. For events 
that contain only a single pp interaction, the location 
of the interaction vertex can be determined from the 
time difference between the hits at the two ends of 
the detector to an accuracy of 3 cm. This interaction 
vertex is used in the last level of the trigger. 

The next trigger level consists of an AND-OR 
decision network programmed to trigger on a pp 
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crossing when several preselected conditions are met. 
This decision is made within the 3.5 [is time inter- 
val between beam bunch crossings. The signals from 
2x2 arrays of calorimeter towers ("trigger towers"), 
covering 0.2x0.2 in rj x </>, are added together elec- 
tronically for the EM sections ( "EM trigger towers" ) 
as well as for all sections, and shaped with a fast rise 
time for use at this trigger level. An additional trig- 
ger processor can be invoked to execute simple al- 
gorithms on the limited information available at the 
time of the AND-OR network. These algorithms use 
the energy deposits in each of the calorimeter trigger 
towers. 

The final software-based level of the trigger con- 
sists of an array of 48 VAXstation 4000 computers. 
At this level, complete event information is available 
and more sophisticated algorithms are used to refine 
the trigger decisions. Events are accepted based on 
certain preprogrammed conditions and are recorded 
for eventual offline reconstruction. 



III. DATA SELECTION 



A. Trigger Filter Requirements 



The trigger requires two electrons which satisfy 
the isolation requirement, each with Et>20 GeV. 
Figure ^| shows the measured detection efficiency of 
the electron filter as a function of p T for a thresh- 
old of 20 GeV. We determine this efficiency using 
Z boson data taken with a lower threshold value 
(16 GeV). (The efficiency corresponds to the frac- 
tion of electrons found at the higher threshold.) The 
curve is a paramete rizatio n used in the simulation 
described in Section [11 D. 
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We require the transverse energy, Et {EsmO), of 
one or more trigger towers to be greater than 10 
GeV. The trigger processor computes an EM trans- 
verse energy by combining the Et of the EM trig- 
ger tower (that exceeded some threshold) with the 
largest signal in the adjacent EM trigger towers, but 
doing this only if the original EM signal has at least 
85% of the energy of the entire trigger tower (includ- 
ing hadronic layers). 

For the accepted trigger tower, a software algo- 
rithm finds the most energetic of four sub-towers, 
and sums the energy in a 3x3 array of calorime- 
ter cells around it. It examines the longitudinal 
shower profile by checking the fraction of the total 
energy found in different EM layers. The transverse 
shower shape is characterized by the pattern of en- 
ergy deposition in the third EM layer. The difference 
between the energies in concentric regions centered 
on the most energetic tower covering 0.25x0.25 and 
0.15x0.15 in r\ x <j> must be consistent with expec- 
tations for an electron shower. The trigger also im- 
poses an isolation condition requiring 



P T 



< 0.15 



(5) 



where the sum runs over all cells within a cone of 
radius 1Z = \/ Acf> 2 + A?/ 2 = 0.4 around the elec- 
tron direction and p T is the transverse momentum of 
the electron, based on its energy and the z-position 
of the interaction vertex as measured by the ho- 
doscopes. 
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FIG. 2. Electron detection efficiency as a function of 
electron Et at the trigger level. The efficiency is essen- 
tially flat above our final Et cutoff of 25 GeV. 



B. Fiducial and Kinematic Requirements 

Events passing the filter requirement are analyzed 
offline where they are reconstructed with finer pre- 
cision. The two highest-i/r electron candidates in 
the event, both having Et> 25 GeV, are used to 
reconstruct the Z boson candidate. One electron 
is required to be in the central region, |?7dct| < 1-1 
(CC), and the second electron may be either in the 
central or in the forward region, 1.5 < |??dct| < 2.5 
(EC). This yields two topologies for the selected 
events: CCCC, where both electrons are detected in 
the central region, and CCEC, where one electron is 
detected in the central region and the other in the 
forward region. In order to avoid areas of reduced 
response between neighboring (f> modules of the cen- 
tral calorimeter, the (f> of any electron is required to 
be at least 0.05 x 27r/32 radians away from the po- 
sition of a module boundary. Finally, the events are 
required to have an invariant mass near the known 
value of the Z boson mass, 75 < M ee < 105 GeV. 
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C. Electron Quality Criteria 



To be acceptable candidates for Z production, 
both electrons are required to be isolated and to sat- 
isfy offline cluster-shape requirements. Additionally, 
at least one of the electrons is required to have a 
spatially matching track associated with the recon- 
structed calorimeter cluster. 

The isolation fraction is defined as 



^1000 



/iso — 



E — E 

Eccivo. 



(6) 



where E cone is the energy in a cone of radius 1Z = 0.4 
around the direction of the electron, summed over 
the entire depth of the calorimeters, and E core is the 
energy in a cone of 1Z = 0.2, summed over only the 
EM calorimeter. Both electrons in the data sample 
are required to have /j SO < 0.15. 

We test how well the shape of any cluster agrees 
with that expected for an electromagnetic shower by 
computing the quality variable Xhm ^ or a ^ ce ^ ener ~ 
gies using a 41-dimensional covariance matrix called 
the H-matrix ]2^] . The covariance matrix is deter- 
mined from GEANT-based simulations 26 2^, which 
were tuned to agree with test beam measurements. 
Both electrons in the sample are required to have a 
tight selection of Xhm < 100- 

The quality of the spatial match between a re- 
constructed track and an electromagnetic cluster is 
defined by the variable 



As 
Is" 



Az 

Sz 



(7) 



where As is the distance between the centroid of the 
cluster in the third EM layer and the extrapolated 
trajectory of the track along the azimuthal direction, 
and Az is the analogous distance in the z-direction. 
For EC electrons, z is replaced by r, the radial dis- 
tance from the center of the detector. The param- 
eters Ss — 0.25 cm, Sz = 2.1 cm, and Sr — 1.0 cm, 
are the resolutions in As, Az, and Ar, respectively. 
At least one of the candidate electrons is required 
to have a < 5 for candidates with \r]det\ < 1-1 and 
a < 10 for candidates with 1.5 < |?7dct| < 2.5. 

The total integrated luminosity of the data sam- 
ple is 111 pb _1 . After applying the selection crite- 
ria, 6407 events remain, with 3594 events containing 
both electrons in the central region and 2813 events 
containing one electron in the central region and one 
in the forward region. Figure [| shows the mass and 
Pt distributions (for 75 < M ee < 105 GeV) in the 
final data sample. There are 157 events with p^>50 
GeV, and the event with the largest pr has pt = 280 
GeV. 
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FIG. 3. (a) Mass distribution for all accepted electron 
pairs, and (b) the pr distribution for those pairs with 
75 < M ee < 105 GeV. 



D. Resolutions and Modeling of the Detector 

Both the acceptance and the resolution-smeared 
theory are calculated using a simulation technique 
originally developed for measuring the mass of the 
W boson J2l| and inclusive cross sections of the W 
and Z bosons psfl , with minor differences arising 
from small differences in the selection criteria. We 
briefly summarize the simulation here. 

The mass of the Z boson is generated according to 
an energy-dependent Breit-Wigner lineshapc. The 
Pt and rapidity (y) are chosen randomly from grids 
created with the computer program legacy |TJ| 
which calculates the Z boson cross section for a 
given pt, y, and Q. For calculating the grids, we 
use a fixed value for the mass of the Z boson of 
91.184 GeV. We match the low-p^ and high-px re- 
gions following the algorithm used in the program 
RESBOS [[l9| to produce a grid of pt and y values, 
weighted by the production cross section, calculated 
to NNLO. The primary vertex distribution for the 
event is modeled as a Gaussian with a width of 27 cm 
and a mean of —0.6 cm, corresponding to the width 
and offset measured in the data. The positions and 
energies of the electrons are smeared according to 
the measured resolutions and corrected for offsets 
in energy scale caused by the underlying event and 
recoil particles emitted into the calorimeter towers. 
Underlying events are modeled using data from ran- 
dom inelastic pp collisions with the same luminosity 
profile as the Z sample. 

The electron energy and angular resolutions are 
tuned to reproduce the observed width of the Z — ► 
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ee mass distribution at the Z resonance. The frac- 
tional energy resolution can be parameterized as 
a function of electron energy as AE/E = C © 
S/\/Et. The sampling term, S, was obtained from 
measurements made in a calibration beam, and is 
0.135 GeV 1 / 2 for the CC and 0.157 GeV 1 / 2 for the 
EC |^,^|. The constant term, C was determined 
specifically for our selection criteria. In the CC, the 
value is C = 0.014 ± 0.002 and in the EC the value is 
C = O.Oig^J. The uncertainty is dominated by the 
statistics of the Z — > ee sample. The uncertainty in 
the polar angle of the electrons is parameterized in 
terms of the uncertainty in the center-of-gravity of 
the track used to determine the polar angle. Fig- 
ure |^ compares electrons from Z boson data with 
simulated results for distributions in electron Et, 
pseudorapidity, and 4>. 

In addition to the smearing of the electron energies 
and positions, certain specific features of the exper- 
iment are also modeled in the simulation in order 
to more closely represent the data. A parameteriza- 
tion of the rise in efficiency of the trigger as a func- 
tion of electron Et is included, as well as a param- 
eterization of the tracking efficiency as a function 
of electron pseudorapidity. Both efficiencies have a 
negligible effect on the shape of the pr distribution. 
Details of the detector simulation can be found in 
Refs. [5TU321. 
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FIG. 4. Comparison of electron Et, r\, and <p, from Z 
boson data (crosses) to results of the detector simulation 
(dashed). 



IV. EFFICIENCY 

We determine the efficiency of the event selec- 
tion criteria as a function of the pt of the Z boson, 
normalizing the result to the integrated total cross 
section for Z boson production as measured at D0 
(a z ■ B(Z -> ee) = 221 pb) @). 

Of all the selection criteria, the electron isolation 
requirement has the largest impact on the observed 
Pt of the Z boson. Nearby jet activity spoils the 
isolation of an electron, causing it to fail the selec- 
tion criteria. The effect depends upon the detailed 
kinematics of the event, in particular, the location 
of hadronic activity (e.g., associated jet production) 
and the pt of the vector boson. 

Two methods have been used to determine the pt- 
dependence of the electron identification efficiency. 
In the first method, the effect of jet activity near 
an electron shower is parameterized in terms of the 
component of the hadronic recoil energy (u) pro- 
jected onto the vector p T . This is denoted as U|| fS3[ . 
The relationship between p\ and u\\ is illustrated in 
Fig. ||. We used a combination of simulated electrons 
and W boson data to obtain the efficiency for identi- 
fying electrons as a function of u\\ . Electron showers 
were generated using the GEANT detector-simulation 
program, and the parameters for the simulated elec- 
trons (e.g., Et, isolation, Xhm) agreed well with 
those observed in W boson data pl| . The agree- 
ment suggests that the effect of hadronic activity on 
the electron is well-modeled in the simulation. Al- 
though our parameterization is obtained using elec- 
trons from W events, we apply it to electrons from 
Z boson events (which have very similar energy dis- 
tributions due to hadronic recoil), because the pa- 
rameterization reflects the effect of hadronic activity 
on high-pT electrons, regardless of the origin of that 
activity. 

The electron identification efficiency as a function 
of U|| is parameterized as: 



a if (uh < u ) 

a(l - s(u|| - u )) if (u\\ > u ) 



(8) 



where u D is the value of tin at which the efficiency 
begins to decrease with un, and s is the rate of de- 
crease. The values obtained from the best fit are 
are u a = 3.85 ± 0.55 GeV and s = 0.013 ± 0.001 
GeV -1 . The parameter a reflects the overall effi- 
ciency, which, as we have indicated, is obtained from 
a normalization to the overall selection efficiency. 
The final event efficiency as a function of pt of the 
Z boson, shown in Fig. ^|, is obtained from the de- 
tector simulation, by comparing the pt distribution 
with and without the it || correction. The final event 
efficiency is insensitive to the use of different param- 
eterizations of the U|| efficiency in the EC versus the 
CC. A more detailed description of the method used 
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FIG. 5. Illustration of the relationship between the 
transverse momentum of the electron, the vector Et of 
the hadron recoil (it) in the calorimeter, and um, the 
projection of the recoil onto the transverse direction of 
the electron. In the particular example illustrated here, 
it|| is negative. 



to obtain the um parameterization can be found in 
Ref. |@. 

In the end, the it|| parameterization of the event 
identification efficiency alone is unsatisfactory for 
application to this measurement. In particular, 
that analysis required < 30 GeV, thereby restrict- 
ing applicability to that region. To obtain a rea- 
sonable parameterization of the electron identifica- 
tion efficiency for all values of pr, we extract the 
Z boson identification efficiency from events gener- 
ated with herwig |34], smeared with the D0 de- 
tector resolutions, and overlaid onto randomly se- 
lected pp collisions ("zero-bias" events). The effi- 
ciency as a function of px is defined by the ratio of 
the pt distribution for events with resolution smear- 
ing and kinematic, fiducial and electron quality re- 
quirements imposed, to that with only kinematic 
and fiducial requirements. Figure ^a compares the 
efficiency as a function of pr using the u\ \ parmeteri- 
zation with that using the detector-smeared herwig 
events. The distributions have been normalized to 
each other in the region px<30 GeV. Figure ^(b) 
shows the ratio of the two normalized results for 
Pt<30 GeV. The agreement of the herwig analysis 
with the Uii analysis is taken as confirmation of the 
validity of the herwig result for all pr- (The model 
for the um analysis has been shown to be reliable for 
p T <30 GeV.) 

In normalizing our efficiency to the previously de- 
termined inclusive Z boson event selection efficiency, 



we use the combined CCCC and CCEC efficiency 
of 0.76 |p8| . We fit the herwig result to a lin- 
ear function in the region pt<18 GeV, and a con- 
stant in the region pt>18 GeV, to obtain the pr- 
dependent event selection efficiency for all pj val- 
ues. The parameterization is shown in Fig. ffl. The 
PT-dependence of the efficiency, in absolute terms, 
is given by e = 0.78 - 0.004p T , for p T < 18 GeV, and 
0.73 for p T >18 GeV. 

We assume that the efficiency above 100 GeV is 
the same as in the region of 18-100 GeV. This is the 
simplest assumption we can make given the statistics 
of the simulation. The efficiency at high-pT cannot 
be greater than at pr = 0, which would correspond 
to about a 1.5 standard deviation change in the cross 
section in that region, and this difference would be 
reflected in the uncertainty on the extracted differen- 
tial cross section. We do not expect the efficiency to 
decrease in the region beyond 100 GeV, because the 
jets in such events will tend to be in the hemisphere 
opposite to the electrons. Events with high jet mul- 
tiplicity may have instances in which the large-i?T 
jets balance most of the transverse momentum of 
the event, but smaller-£/r jets can overlap with one 
of the electrons. However, because the electrons are 
very energetic, low energy jets are not likely to af- 
fect the efficiency of the isolation criteria. We assign 
estimated uncertainties on the efficiency of ±3% in 
the bin below 18 GeV, and ±5% in the region above 
18 GeV. 
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FIG. 6. (a) Comparison of the Z boson selection effi- 
ciency as a function of pr as determined using herwig 
(dashed crosses), and as determined using a parameter- 
ization of the single-electron efficiency as a function of 
um (solid crosses), (b) The ratio of the two methods in 
the range 0-30 GeV, where they are expected to agree. 
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FIG. 7. Final event identification efficiency as a func- 
tion of pr, based on herwig events; the line is the pa- 
rameterization used in calculating the final cross section. 



V. ACCEPTANCE 

The parameterized detector simulation referred to 
in Section [II is used to determine the overall accep- 
tance as a function of pt of the Z boson. The effects 
of the trigger turn-on in Et, the rapidity cut-offs, 
the <f> module boundaries in the central calorime- 
ter, the pseudorapidity dependence of the tracking 
efficiency, and the final Et requirements are all in- 
cluded in the calculation of the acceptance. Figure || 
shows the relative effects of the requirements on the 
electron Et and pseudorapidity, and of the trigger 
and tracking efficiency on the acceptance as a func- 
tion of pt- As can be seen, the strongest effects 
come from the electron Et and pseudorapidity re- 
quirements. The dip in relative acceptance seen in 
Fig. |8|(a) for middle values of pt results from one 
of the electrons carrying most of the pt of the Z 
boson-one electron can have a relatively large Et 
while the other has relatively small Et- However, 
as the pt of the Z boson increases beyond 45 GeV, 
this asymmetry is no longer allowed - both electrons 
must have relatively large Et- The monotonic rise 
of the relative acceptance in Fig. [|(b) is due to the 
increasing "centrality" of the event-as pt increases, 
the rapidity of the Z boson is closer to zero. As can 
be seen in Figs. ||(c) and (d), the imposition of the 
other selection criteria merely changes the normal- 
ization and does not affect the shape as a function 
of pt- 

The mass requirement on the dielectron pairs has 
been ignored in the final acceptance calculation. 
Figure Q compares the pt distribution for dielectron 



pairs with invariant mass near that of the Z boson 
to those with invariant mass above and below the 
nominal Z boson mass, and supports the expecta- 
tion that any pt dependence on mass (near the Z 
boson mass peak) is very small. 

Figure [To] shows the acceptance for the CCCC and 
CCEC event topologies, as well as for the combined 
event sample. Here we see the increased centrality of 
the events as a function of pt, noting the increasing 
acceptance for the CCCC events in contrast to the 
decreasing acceptance for the CCEC events. The dip 
and rise in Fig. 10|(b) are due to competing effects of 
the electron Et and pseudorapidity requirements. 

The effect of uncertainties in the energy scale and 
resolution, the tracking resolution, and the trigger 
efficiency is assessed for each bin of pt by vary- 
ing the values of these parameters by their mea- 
sured uncertainties. Figure ^ shows the nominal 
acceptance and those obtained by varying the val- 
ues of the parameters. The largest differences are 
observed at high pt- If we parameterize this sys- 
tematic uncertainty as a linear function of pt, we 
obtain S acc = ±(0.01 + 0.0001p T )- This resulting 
band of uncertainty is also shown in Fig. [ll]. 

Because we determine the acceptance bin by bin 
in pt, we are relatively insensitive to the underly- 
ing model for the pt spectrum used in the detector 
simulation. Nevertheless, we are sensitive to the as- 
sumed rapidity distribution of the Z boson in each 
bin of pt- The uncertainty in the predicted rapidity 
of the Z boson is expected to be dominated by the 
uncertainty in the pdf's used for modeling Z pro- 
duction. The uncertainty in acceptance due to the 
choice of pdf has been found to be ±0.3% for the 
inclusive measurement of the Z boson cross section 
p8[ . This constrains the uncertainty in the low-pT 
region, where the cross section is largest, to a value 
that is far smaller than the uncertainty from varia- 
tions in the parameters of the model of the detec- 
Figure 



tor 
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shows that the rapidity distributions 
obtained from the detector simulation and for data 
agree for both low and high values of pt', we there- 
fore ignore any additional uncertainty in the accep- 
tance due to the modeling of the rapidity of the Z 
boson. 



VI. BACKGROUNDS 

The primary background to dielectron production 
at the Tevatron is from multiple-jet production from 
QCD processes in which the jets have a large elec- 
tromagnetic component (most of the energy is de- 
posited in the EM section of the calorimeter) or they 
are mismeasured in some way that causes them to 
pass the electron selection criteria. There are also 
contributions to the Z boson dielectron signal that 
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FIG. 8. Effect of requirements on (a) Et, (b) r), (c) 
trigger Et cut-off, and (d) tracking efficiency on the rel- 
ative acceptance as a function of pr- 
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FIG. 10. Final acceptance as a function of pr for (a) 
the CCCC, (b) CCEC event topologies, and (c) for the 
combined event sample. 
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FIG. 9. Comparison of the transverse momentum dis- 
tribution of dielectron pairs with mass very close to the 
nominal Z boson mass, 90 < M ee < 92 GeV (solid cir- 
cles) to those in the mass regions 75 < M ee < 90 GeV 
(open circles) and 92 < M ee < 105 GeV (closed trian- 
gles). 
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tion of the uncertainty as a function of pr, as discussed 
in the text. 
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(a) all values of p T , (b) p T <20 GeV, and (c) 20<p T <50 
GeV. 



are not from misidentification of electrons, but cor- 
respond to other processes that differ from the one 
we are trying to measure, e.g., Z — > t + t~ and tt 
production. Such processes are irreducible due to 
the fact that they have the same final event signa- 
ture as the signal, and often have pr dependences 
that can differ from the Z/"f* mediated production 
of the Z boson and Drell-Yan pairs. These must be 
determined and accounted for in any comparison of 
data with theory. 

Both the normalization and the shape of the mul- 
tijet background as a function of px are determined 
from data. Three types of backgrounds have been 
studied to examine whether differences in produc- 
tion mechanism or detector resolution would pro- 
duce a significant variation in the background: di- 
jet events (from multijet triggers), direct-7 events 
(from single photon triggers), and dielectron events 
in which both electrons failed the quality criteria 
(from the Z boson trigger). For the dijet events, we 
selected the two highest- Et jets and reconstructed 
the "Z boson" as if the jets were electrons. Simi- 
larly, for the direct-7 events, we selected the highest- 
Et photon candidate and the highest-£r jet in the 
event. For the failed-dielectron sample, we used the 
two highest-£/r electron candidates whose cluster 
shape variable (Xhm) did not match well with that 
of an electron. For all three backgrounds, the "elec- 
tron" objects were required to satisfy the same Ex 
and 77 criteria as the data sample. 

Figures p~3| — p~6| show the invariant mass and pt dis- 
tributions for the background samples in both the 



CCCC and CCEC event topologies. The direct- 
photon and failed-dielectron events agree in the mass 
and pt distributions. The Kolmogorov-Smirnov 
probability (Pks) for the two mass distributions is 
0.78 and between the two px distributions it is 0.97. 
Thepr distribution from the dijet sample also agrees 
well with the direct-7 an d failed-dielectron samples, 
with Pks — 0.51 and Pks = 0.57, respectively. The 
dijet mass distribution does not agree as well, giv- 
ing Pks = 0.005 when comparing to the direct-7 
sample and Pks = 0.1 when comparing to the failed- 
dielectron sample. The difference is likely due to the 
poorer jet-energy resolution compared to the elec- 
tron energy resolution. This difference in the shape 
of the invariant mass is included in the systematic 
uncertainty on the backgro und n ormalization, and 
is a small effect (see Section VTA ) . 
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FIG. 13. Invariant mass distributions for the three 
types of QCD multijet background samples in the CCCC 
topology: (a) dijet data sample, (b) direct-7 data sam- 
ple, (c) failed-dielectron data sample. We show the dis- 
tributions for all three data samples in (d). 



A. Multijet Background Level 

Because the mass distribution for the multijet 
background samples depends on event topology, the 
level of the multijet background is determined sepa- 
rately for CCCC and CCEC dielectron events. Us- 
ing this background, and the contribution from the 
Z boson, we can obtain the relative background 
fraction through a maximum-likelihood fit for the 
amount of background and signal in the data. 

We use the PYTHIA event generator (3^] to produce 
the invariant mass spectrum for the signal. Con- 
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FIG. 14. Invariant mass distributions for the three 
types of QCD multijet background samples in the CCEC 
topology: (a) dijet data sample, (b) direct-7 data sam- 
ple, (c) failed-dielectron data sample. We show the dis- 
tributions for all three data samples in (d). 
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FIG. 15. Comparison of shapes of transverse momen- 
tum distributions for the three multijet background sam- 
ples for the CCCC event topology. 
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FIG. 16. Comparison of shapes of transverse momen- 
tum distributions for the three multijet background sam- 
ples for the CCEC event topology. 

tributions from both Z boson and Drell-Yan pro- 
duction and their quantum-mechanical interference 
are included in the calculation. The generated four- 
momenta are smeared using the detector simula- 
tion described previously. We obtain the amount 
of multijet background in the data by performing 
a binned maximum-likelihood fit to the sum of the 
signal (pythia) and background: 

^Vdata(TOz) = CiNp YTmA (m.i) + C 2 ^backgrounding) 

(9) 

where c\ and ci are the normalization factors for 
the signal and background contributions, respec- 
tively, and rrii is the zth mass bin. The fit was per- 
formed in the dielectron invariant mass window of 
60 < M ee < 120 GeV. Figure [l?] shows the best fit to 
the dielectron invariant mass, separately for CCCC 
and CCEC topologies using the direct-7 sample as 
the background. Using the other two background 
samples yields similar results. The final value for the 
fraction of multijet background in the data, /back, is 
defined by normalizing the fit parameter c\ to the 
number of events observed in the mass window of 
the Z boson (75 < M ee < 105 GeV): 



/1 



back 



where 



Cl 



iVtotai(data) iV 75 _ 105 (background) 



A/totai (background) iV 75 _i 5(data) 



(10) 



-Ntotal (sample) = ^ ^sample (w,) (11) 



all 
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^75- 105 (sample) = ^ N aamp \ e (mi). 

75<mi<105 (12) 

We use the direct-7 sample for the central value 
of the level of multijet background, and use the sta- 
tistical uncertainty from that fit. We also assign 
a systematic uncertainty associated with our choice 
of mass window used in the fit and for differences 
in the background models. We assign a systematic 
uncertainty to the background normalization that 
corresponds to half of the maximum difference from 
the central value in the determined background frac- 
tions. The background values for each topology and 
the resulting uncertainties are summarized in Ta- 
ble § 

Combining the uncertainties in quadrature, we ob- 
tain a background fraction of (2.45 ± 0.50)% for the 
CCCC topology and (7.09 ± 1.00)% for the CCEC 
topology. Weighting the background fractions by the 
relative number of events in each topology, we obtain 
a total multijet background level of (4.44±0.89)%. 




110 120 130 
M(GeV) 



110 120 130 
M(GeV) 

FIG. 17. Comparison of the dielectron invariant mass 
distribution (closed circles) and the background (dashed 
line) to the fit to pythia Z/"/* and background (solid 
line) for (a) CCCC and (b) CCEC topologies. 



B. pT-Dependence of the Multijet Background 

The direct-7 sample is used to determine the 
shape of the background distribution for several rea- 
sons. First, this sample has the greatest number of 
events. Second, we expect the direct-7 data sam- 
ple to provide a good approximation of the combi- 
nation of backgrounds from dijet and true direct-7 
production because about half of the direct-7 sam- 
ple consists of misidentified dijets, and therefore has 



Background Model 


CCCC 


CCEC 


direct-7 
55 < M ee < 125 
65 < M ee < 115 

dijets 
failed dielectrons 


(2.45 ± 0.41%) 
(2.10 ±0.36%) 
(2.74 ±0.51%) 
(1.98 ± 0.35%) 
(2.10 ±0.37%) 


(7.09 ± 0.87%) 
(7.52 ± 0.83%) 
(6.84 ± 0.96%) 
(6.37 ± 0.80%) 
(6.22 ± 0.78%) 


model uncertainty 
window uncertainty 


0.24% 
0.17% 


0.44% 
0.22% 



TABLE I. Background fractions in the two primary 
event topologies. The values in the first five rows include 
only statistical uncertainties for each method. The sys- 
tematic uncertainties obtained by considering variations 
in the background selection and fitting criteria are shown 
in the last two rows. 



Parameter 


p T < 8 GeV 


Pt > 8 GeV 


a 
a 
b 


-0.31 ± 0.02 
0.001 ± 0.002 
0.65 ±0.18 


-0.1014 ± 0.0015 
0.08 ± 0.03 
0.136 ±0.004 



TABLE II. Values of the parameters obtained from 
the fit to the direct-7 background. 



the approximate balance of dijet and direct-7 events 
expected from QCD sources. Third, since events in 
the direct-7 often contain at least one good electron- 
like object, detailed differences between choosing 
electron-like objects and jet objects for reconstruct- 
ing the "Z" boson are smaller here. 

The final shape of the background is obtained by 
combining the CCCC and CCEC samples, weighted 
by the relative contributions to the background. To 
facilitate later analysis, the shape is parameterized 
as a function of px using the following functional 
forms: 



a{ PT + b) 2 e apT if (p T < 8 GeV) 
a{±) 2 + be apT if \p T > 8 GeV) 



(13) 



The function is normalized to be a probability dis- 
tribution, that is the product of the function and the 
total number of background events results in the dif- 
ferential background in each bin of px- Figure [l8| 
shows the results of the fit to the background and 
Table O shows the values of the fit parameters. 



C. Other Sources of Dielectron Signal 

Although Z — > ee and QCD multijet events make 
up nearly all of observed dielectron signal, there 



15 




(a) 



^ 10 

x - 2 

& 10 
§10 

5 10 

"S - 5 
g 70 
a, -6 
10 

-7 

10 




8 10 
PjtGeV) 

(b) 







50 



100 



150 



200 



250 300 
PjfGeV) 

FIG. 18. The results of the fit to the background as a 
function of p T . For (a) p T <8 GeV, the xVd.o.f. = 0.9, 
and (b) for p T >8 GeV, the x 2 / d -°- f - = 0.7. 



are contributions from other sources, such as Z —* 
t+t~, it, and diboson {WW, ZZ, WZ, W~f, Z~f) 
production in dielectron final states. The expected 
contributions from these sources are estimated be- 
low. 

The dielectron event rate from Z — > t + t~ produc- 
tion in our accepted mass range is calculated to be 
< 2.6 x 1CP 6 per Z — > ee event |3(|. The events were 
generated with the herwig simulator and smeared 
with the D0 detector resolutions. For the current 
sample, this corresponds to less than 0.009 events 
for all values of dielectron pp. We therefore ignore 
this contribution to the signal. 

The dielectron background contribution from it 
production is concentrated at high p^- The frac- 
tion was determined using the herwig simulator for 
it production, smeared with the known D0 detec- 
tor resolutions. Electron contributions from both 
W — > ev and W — + tX — + eX channels were con- 
sidered. For a it cross section of 6.4 pb [^Tj, and 
the standard branching ratios for the W boson, the 
calculated geometric and kinematic acceptance from 
herwig is 0.01±0.006. Including electron identifica- 
tion efficiency for dielectron events, we expect about 
0.36 events in the entire sample and about 0.2 events 
with dielectron pt> 50 GeV. Considering the small 
number of events expected, the it contribution is also 
ignored. 

We considered WW, ZZ, WZ, and Wy events 
generated with the herwig simulator and smeared 
with the known D0 detector resolutions. All of these 
backgrounds are small, and we therefore focus on 
any possible effects on our measurement at high-pr, 



Process 


Acceptance 


a ■ B (pb) 


Expected 
events 


tt 


0.01±0.006 


0.43±0.01 


^0.2 


Wy 


<0.0003 


11.3±0.3 


< .3 


WW 


0.016±0.007 


0.12±0.03 


wO.15 


WZ 


0.016±0.007 


0.08±0.01 


wO.l 


ZZ 


0.046±0.002 


0.03±0.01 


^0.005 



TABLE III. The expected number of events from di- 
boson and top quark processes. Each channel assumes a 
total luminosity of 108.5 pb -1 , and a dielectron identifi- 
cation efficiency of 0.73. 



where there are relatively few events and effects of 
even a small background contamination could be sig- 
nificant. 

The resulting acceptances and expected number 
of background events with px>50 GeV are given in 



Tabic 



III 



No W*y events out of approximately 3000 
generated passed the selection requirements, because 
very few such events have photons with Et>25 GeV 
and an invariant mass (M e7 ) near the Z boson mass. 
The table includes the assumed production cross 
sections multiplied by branching ratios (a ■ B) for 
W and Z boson into electron states. The WW 
cross section (10. 21^ pb) and branching ratio to 



|. The 
is ob- 



dielectrons (0.011) are obtained from Ref. |3 
value a(Wj) ■ B(W -> ev) = 11.3±^ pb 
tained from Ref. |3^| , and assumes p7p > 10 GeV 
0.7. The standard model WZ cross 



and Ai? e7 > 



and the ZZ cross 
Given their small 



section is taken from Ref. 3£j 
section is taken from Ref. [40 
size, all of these contributions have been ignored in 
our analysis. 



VII. MEASURED da/dp T 



Table IV shows the values for each of the individ- 
ual components of the measurement: the number of 
events observed for each bin of pt, the product of 
the efficiency and the acceptance (e • a) , and the ex- 
pected number of background events (b) . The associ- 
ated uncertainties are also included. We combine the 
geometric acceptance and identification efficiencies 
into a single overall event efficiency by taking their 
product. We assume that the uncertainties are well- 
described as Gaussian distributions and add them in 
quadrature to obtain the uncertainty S(e ■ a). 

The measured differential cross section, da'/dpr, 
is obtained by calculating the cross section in each 
bin of pt, accounting for the effects of efficiency, 
acceptance, and background, but not accounting for 
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the effects of detector smearing. That is, 



(14) 



where cr^ is the measured cross section in bin i and 
Af m is the width of the bin in pt- 

We obtain the cross section and uncertainty in 
each bin using the methods of statistical inference. 
We relate the expected number of events y, in each 
bin to the underlying cross section JIT 42 1: 



jj, = £ecr' + b, 



(15) 



where £ is the total integrated luminosity, e is the 
overall detection efficiency for the process, and b is 
the number of background events. A value of /i is 
determined for each bin of px- 

We relate the observed number of events and the 
expected number of events through a probability dis- 
tribution, in our case an assumed Poisson distribu- 
tion, 



P(d\a',e,b,£,I) 



e-( Cea ' +b \£e<j' + b) d 



(16) 



where d is the number of events observed and / refers 
to the assumptions implicit in deriving the probabil- 
ity density |42| . 

Applying B ayes' Theorem, we invert the probabil- 
ity in Eq. fLq, 



P{a',e,b,£\d,I) = 

P(d\d',e,b,C,I)P(<T',e,b,C\I) 

c 



(17) 



where £ normalizes the probability such that 
/ P(a' ,e,b,£\d,I)dfl = 1, where dfl denotes 
that the integration is over all relevant variables. 
P(<r' , e, b, £\d, I) is the joint posterior probability, 
describing the probability of a particular set of 
cr',e, b, and £ given the results from our data. 
P(er', e, b, £\I) is the joint prior probability, describ- 
ing the probability of a particular set of a', e, 6, £ be- 
fore taking our data into account. P(d\a' ', e, 6, £, I) 
is the likelihood function for our data. Assuming 
that the individual parameters arc logically indepen- 
dent, e.g., the cross section does not depend on the 
background, then Eq. 17 can be rewritten as 



P{a',e,b,£\d,I) = 
P(d\a>, e, b, £, I)P{a>\I)P{e\I)P{b\I)P(£\I) 

c 



(18) 



We are not interested in the values of the param- 
eters e, b, and £, and we eliminate the dependence 
of the joint posterior probability on these nuisance 
parameters by integrating over their allowed values, 
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FIG. 19. The measured differential cross section, not 
corrected for detector smearing, (a) for pt<50 GeV and 
(b) for all pt- 



a process called marginalization. To extract our 
results, we calculate 



P(a'\d,I) = J dbd£deda' 



P(d\a\ e, b, C, I)P(a'\I)P(e\I)P(b\I)P(£\I) 

c 



(19) 



In the calculation of the binned differential cross sec- 
tion, the uncertainty on the integrated luminosity 
changes only the overall normalization of the distri- 
bution, which is already accounted for in our nor- 
malization to the D0 measurement of &z^ee- We 
therefore use a delta function as the prior for the 
integrated luminosity distribution. We assume the 
priors for the efficiency (e) and background (b) to be 
Gaussian distributed, with their estimated mean val- 
ues and standard deviations as the means and widths 
of the Gaussians. The prior probability distribution 
for the cross section in each bin is taken to be in- 
dependent of er (uniform for the range [u m in,^max] 
where cr m j n > 0) and the total range is at least ±6 
standard deviations around the mean. 

The integration in Eq. 19 is performed using the 
numerical integrator MISER [j43| , and the results are 
shown in Fig. [20|. Since the probability distributions 
for all but the highest-p^ bin are nearly Gaussian, we 
assign the final value of the cross section for each bin 
in pt to be the mean of the probability distribution 
with uncertainties set equal to the standard devia- 
tion about the mean. For the highest-py bin, we use 
the most probable value for the cross section with 
upper and lower uncertainty values circumscribing 
the narrowest 68% confidence interval. The integral 
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over px of the differential distribution is normalized 
to the inclusive cross section for Z boson production 
measured by D0 |§. Table [fv| gives the values of 
the measured differential cross section in each bin 
of pt, not corrected for detector smearing, and Fig- 
ure |l^ displays the results as a function of pr- 



VIII. FIT TO NON-PERTURBATIVE 
PARAMETERS 

As discussed in Section |, the current theoretical 
understanding of the pr distribution of Z bosons 
uses fixed-order perturbative calcuations to describe 
the high-pT region and resummation calculations 
of the perturbative solution to describe the low- 
Pt region. At the smallest values of pr, a pa- 
rameterization must be invoked to account for non- 
perturbative effects that are not calculable in per- 
turbative QCD. The generic form for the function is 
given in Eq. ^, however, one must choose particular 
functional forms for hi(x,b) and h 2 (b). Historically 
there are two versions for the choice of this parame- 
terization. The first, from Davies, Weber, and Stir- 
ling (DWS) @, has the form: 



S 



DWS 
NP 



{b 1 Q 2 )^g 1 b 2 +g 2 b 2 \n^-). (20) 



The values of g± and g 2 are determined by fitting to 
low-energy Drell-Yan data, yielding gi — 0.15 GeV 2 
and g 2 = 0.4 GeV 2 , where Q = 2 GeV and b max — 
0.5 GeV" 1 (see Eq. |). They used the pdf 's of Duke 
and Owens fig] . The second is from Ladinsky and 
Yuan [HI: 



makes a far larger contribution to the value of S^p. 
The relative importance of g 2 over g\ comes from 
the ln(Q 2 /Q 2 ) term. 

Because the width of the Z boson is « 2.5 GeV, for 
purely phenomenological needs the non-perturbative 
physics can be parameterized using a single parame- 
ter g' = gi + g 2 log(M|/Q 2 ) Q. However, because 
the general form of Snp is theoretically motivated, 
we preserve the form of Eq. [2l], focusing on the value 
of g 2 , the parameter we are most sensitive to. 

We perform a minimum-x 2 fit to determine the 
best value of g 2 from our data. For the purposes of 
the fit, we fix g x = 0.11 GeV 2 and g 3 = -1.5 GeV" 1 , 
as suggested by Ladinsky and Yuan |ll| . We use the 
program LEGACY Q with the CTEQ4M pdf 's @ 
to generate the da' /dpr distribution for the Z bo- 
son and match the low-p^ and high-py regions using 
the prescription in RESBOS, obtaining a single grid 
for all values of pr calculated to NNLO. We smear 
the prediction with the D0 detector resolutions and 
fit the resulting px distribution to our measured re- 
sult. The x 2 distribution as a function of g 2 is well- 
behaved and parabolic and when fit to a quadratic 
function yields a value of 0.59±0.06 GeV 2 at the 
minimum, with x 2 /d.o.f= 10.6/10. 

For completeness, we also fit for the individual 
values of g\ and g 3l using the Ladinsky and Yuan 
values for the two parameters not being fitted. The 
variation in \ 2 of gi and g 3 are also well-behaved 
and parabolic, and the fit yields gi = 0.09 ± 0.03 
GeV 2 and g 3 = -1.1 ± 0.6 GeV" 1 . The value of 
gi agrees with the Ladinsky- Yuan result, and is of 
comparable precision. The value of g 3 also agrees 
with the Ladinsky- Yuan result, but is far less precise. 



S^ Y P (b, Q 2 ) = 

9l b 2 + g 2 b 2 \n(^) + g 1 g 3 bln{WQx i x j ) (21) 

where Xi and the momentum fractions of the 

colliding partons. The values of g\, g 2 , and g 3 are 
determined by fitting to low-energy Drell-Yan data 
and a small sample of Z — > ee data from the 1988-89 
run at CDF (l|, yielding g x = 0.11±g;g| GeV 2 , g 2 = 
Q-58±gi GeV 2 , and g 3 = -1.5±g;J GeV~\ where 
b m ax = 0.5 GeV" 1 and Q a = 1.6 GeV. They used 
the CTEQ2M pdf 's (46) in the fits. 

The Z boson px distribution is by far most sen- 
sitive to the value of g 2 . For measurements at the 
Tevatron at Q 2 = M§, the calculation is nearly in- 
sensitive to the value of g 3 , and only slightly sensi- 
tive to the value of g\. For 173, this insensitivity is 
due to the high energy of the pp beam relative to 
the Q 2 being probed. For a center-of-mass energy of 
s = XiXjS and s « M%, we see that for a measure- 
ment at the Tevatron (y/s — 1.8 TeV), the g 3 term 
becomes <?i(73&ln(100,s/s) « —lAgig 3 b. The g 2 term 
varies as g 2 b 2 \n(Mz /Qo) ~ 4.7c/2& 2 , and therefore 
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Bin 


Pt range 


number 






b 


6b 


da' /dpT 


8{da'/dp T ) 




da/dpr 


S(da /dpr) 


number 


(GeV) 


of events 


e • a 


6{e ■ a) 


(events) 


(events) 


(nb/GeV) 


(nb/GeV) 


a(p T ) 


(nb/GeV) 


(nb/GeV) 


1 


0-1 


156 


0.351 


0.011 


3.28 


0.7 


5.10 


0.45 


1.185 


6.04 


0.53 


2 


1-2 


424 


0.347 


0.011 


8.14 


1.6 


14.0 


0.82 


1.160 


16.2 


0.96 


3 


2-3 


559 


0.346 


0.011 


12.7 


2.5 


18.4 


0.99 


1.108 


20.4 


1.1 


4 


3-4 


572 


0.343 


0.011 


16.1 


3.2 


18.9 


1.0 


1.042 


19.7 


1.1 


5 


4-5 


501 


0.343 


0.011 


18.0 


3.6 


16.4 


0.93 


0.988 


16.2 


0.92 


6 


5-6 


473 


0.342 


0.011 


18.8 


3.8 


15.5 


0.90 


0.965 


15.0 


0.87 


7 


6-7 


440 


0.336 


0.011 


18.5 


3.7 


14.6 


0.88 


0.960 


14.1 


0.84 


8 


7-8 


346 


0.335 


0.011 


17.5 


3.5 


11.5 


0.76 


0.967 


11.1 


0.73 


9 


8-9 


312 


0.334 


0.011 


16.9 


3.4 


10.3 


0.71 


0.972 


10.0 


0.69 


10 


9-10 


285 


0.330 


0.011 


15.2 


3.1 


9.55 


0.69 


0.972 


9.29 


0.67 


11 


10-12 


439 


0.324 


0.017 


26.1 


5.2 


7.46 


0.56 


0.972 


7.25 


0.54 


12 


12-14 


326 


0.317 


0.017 


21.3 


4.3 


5.63 


0.46 


0.967 


5.45 


0.44 


13 


14-16 


258 


0.306 


0.017 


17.4 


3.5 


4.61 


0.41 


0.964 


4.45 


0.39 


14 


16-18 


203 


0.302 


0.016 


14.2 


2.8 


3.67 


0.34 


0.963 


3.54 


0.33 


15 


18-20 


181 


0.297 


0.016 


11.6 


2.3 


3.35 


0.32 


0.958 


3.21 


0.31 


16 


20-25 


287 


0.289 


0.016 


20.5 


4.1 


2.16 


0.19 


0.954 


2.06 


0.18 


17 


25-30 


174 


0.278 


0.015 


12.3 


2.5 


1.37 


0.14 


0.945 


1.29 


0.13 


18 


30-35 


124 


0.270 


0.016 


7.46 


1.5 


1.02 


0.12 


0.944 


0.962 


0.11 


19 


35-40 


104 


0.263 


0.014 


4.51 


0.90 


0.892 


0.10 


0.941 


0.840 


0.10 


20 


40-50 


92 


0.264 


0.014 


4.38 


0.88 


0.392 


0.048 


0.952 


0.373 


0.045 


21 


50-60 


61 


0.274 


0.015 


1.63 


0.33 


0.258 


0.037 


0.974 


0.251 


0.036 


22 


60-70 


40 


0.283 


0.016 


0.616 


0.12 


0.167 


0.028 


0.975 


0.163 


0.027 


23 


70-85 


20 


0.300 


0.017 


0.308 


0.062 


0.054 


0.016 


0.989 


0.053 


0.012 


24 


85-100 


13 


0.319 


0.018 


0.095 


0.019 


0.034 


0.010 


0.988 


0.034 


0.009 


25 


100-200 


15 


0.366 


0.022 


0.130 


0.026 


0.0051 


0.0013 


0.994 


0.0050 


0.0013 


26 


200-300 


2 


0.530 


0.034 


0.038 


0.0076 


0.0004 


+0.00038 
-0.00029 


0.994 


0.0004 


+0.00038 
-0.00029 



TABLE IV. Summary of the results of the measurement of the pr distribution of the Z boson. The range of pr 
corresponds to the intervals used for binning the data. The product of the acceptance and efficiency is given as e • a, 
b is the estimated number of background events, the measured differential cross section is da' /dpr, the correction 
for resolution-smearing is specified by a(pr), and the corrected differential cross section is specified by da /dpr- The 
uncertainty on the differential cross section includes both systematic and statistical uncertatinties, but does not 
include overall normalization uncertainty due to the luminosity of ±4.4%. 
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FIG. 20. Normalized probability distributions for cross sections (<r) in individual bins of pr- The vertical scale in 
each case is Probability /pb. 
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IX. SMEARING CORRECTIONS 

The results shown in Fig. |l9] still contain the resid- 
ual effects of detector smearing. We correct the 
measured cross section for the effects of detector 
smearing using the ratio of generated to resolution- 
smeared ansatz px distributions: 



a{p T ) 



F{pt]92) 



Jdp T R(p T ,p' T )F(p T ;g 2 ) 



(22) 



where p' T is the smeared value of px, ot(px) is the 
correction factor, F(px',g 2 ) is the ansatz function 
with parameter g 2 and R{px,p' T ) is the resolution 
function. 

As the ansatz function, we use the calculation 
from legacy fixing gi — 0.11 GeV 2 and g 3 = -1.5 
GeV -1 . We use g 2 = 0.59 GeV 2 for our central 
value. 

Figure |l] shows the smearing correction as a func- 
tion of pr- The largest effect occurs at low-pr where 
the smearing causes the largest fractional change in 
Pt and where the kinematic boundary at px— re- 
sults in non-Gaussian smearing — the px is preferen- 
tially increased by the smearing rather than being a 
symmetric effect. Table IV includes the value of the 



smearing correction for each bin of px ■ 

It is important that the smearing correction be 
insensitive to significant variations in the ansatz 
function used to generate the correction. We ex- 
amine this issue by varying the parameter g 2 = 
0.59 ± 0.06GeV 2 in the ansatz function by ±1ct, ob- 
taining a variation of <1% for all values of px- For 
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FIG. 22. Fractional uncertainty in the cross section as 
a function of pr including the statistical and systematic 
uncertainties (closed circles) and including only statisti- 
cal uncertainties (open circles) for (a) pr < 50 GeV and 
(b) p T > 50 GeV. 



this variation in the parameter, the ansatz function 
varies by ssl0%. It is useful to compare the level 
of uncertainty in the smearing correction to other 
components of uncertainty in the measurement. Fig- 
ure |22| shows the fractional uncertainty on the dif- 
ferential cross section as a function of px- Both the 
total uncertainties (in which systematic uncertain- 
ties on the background, efficiency, and acceptance 
are included with the statistical uncertainty), and 
the statistical uncertainties alone are shown. The 
variations in the smearing correction are at least a 
factor of five smaller than the other uncertainties 
and therefore can be ignored. 

The uncertainty in the smearing correction is also 
affected by the uncertainty on the values of the res- 
olutions used to generate the smearing. We examine 
this uncertainty by varying the detector resolutions 
by ±1 standard deviation from the nominal values. 
Again, the effect on the smearing correction is negli- 
gible relative to the other uncertainties in the mea- 
surement and this source of uncertainty has been 
ignored. 



X. RESULTS 

Table |^ shows the final numerical results for the 
measurement of da/dpx using a total of 6407 events. 
The uncertainties on the data points include statis- 
tical and systematic contributions. There is an addi- 
tional normalization uncertainty of ±4.4% from the 
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FIG. 23. Plot of the differential cross section (circles) 
as a function of pr- The line is the result of the NNLO 
calculation from RESBOS using the CTEQ4M pdf. The- 
ory has been normalized to the data. 
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FIG. 24. Fractional difference between the data and 
the NNLO calculation from RESBOS for the differential 
cross section as a function of pr- 



uncertainty on the integrated luminosity |2q| that 
is included in neither the plots nor the table, but 
must be taken into account in any fits requiring an 
absolute cross section. 

Figures ^3] ^B] and |2?] show the final, smearing- 
corrected px distribution compared to the fully 
resummed calculation as calculated by RESBOS. 
The calculation uses the CTEQ4M pdf's and 
the Ladinsky-Yuan parameterization for the non- 
perturbative function with the published values for 
the Qi parameters: g\ =0.11 GeV 2 , g 2 = 0.58 GeV 2 
and g 3 = —1.5 GeV -1 . The data points in Figs. g3| 
and |25| are placed at the average px of the bin, as 
given by theory, rather than at the center of the bin. 
(Only the first and twenty-fifth bins change appre- 
ciably, —24% and —17% respectively, relative to the 
bin center.) For the resummed calculation, the total 
cross section predicted by the theory (220 pb) has 
been normalized to the data. A \ 2 comparison of the 
prediction to the data yields a x 2 /d.o.f.= 26.7/25. 
We observe good agreement with the fully resummed 
calculation for all values of px ■ 

Figures |23 and |2^ compare the data to the fixed- 
order (NLOj perturbative calculation as calculated 
by RESBOS and using the CTEQ4M pdf's. We ob- 
serve strong disagreement at low-px, as expected due 
to the divergence of the NLO calculation at px= 0, 
and a significant enhancement of the cross section 
relative to the prediction at moderate values of pr, 
confirming the increase in the cross section from soft 
gluon emission. 



XI. CONCLUSIONS 

We have measured the differential cross section as 
a function of the transverse momentum of the Z bo- 
son. Fitting for the value of the non-perturbative 
parameter g 2 , we obtain g 2 = 0.59 ± 0.06 GeV 2 , 
which is significantly more precise than previous de- 
terminations. We observe good agreement between 
the measurement and the resummation calculations 
for all values of px ■ 
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lation from RESBOS using CTEQ4M as the pdf. Theory 
has been normalized to the data. 
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